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Abstract 

We derive an algorithm for automatic calculation of perturbative /3-functions and anoma- 
lous dimensions in any local quantum field theory with canonical kinetic terms. The infrared 
rearrangement is performed by introducing a common mass parameter in all the propagator 
denominators. We provide a set of explicit formulae for all the necessary scalar integrals up 
to three loops. 
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1. Introduction 

Renormalization group equations are a fundamental tool in modern quantum field theory. 
In phenomenological applications, their evaluation with sufficient accuracy often requires 
finding multiloop contributions to /3-functions and anomalous dimensions. In the present 
paper, we describe a simple algorithm for calculating these quantities in the framework of 
dimensional regular ization and the MS (or MS) scheme. 

In a mass-independent renormalization scheme, /3-functions and anomalous dimensions 
are simply related to coefficients at counterterms which renormalize ultraviolet divergences. 
A remarkable feature of the MS'-scheme is the fact that in its framework all the UV counter- 
terms are polynomial both in momenta and in masses [0J .0 Consequently, a certain expansion 
in external momenta and masses can be performed before integration over loop momenta, 
which radically simplifies the integrals one needs to calculate. 

The main difficulty in this procedure is appearance of spurious infrared divergences. The 
classical method of avoiding them is called "infrared rearrangement" [0, |3|. It amounts to 
adding artificial masses or external momenta in certain lines of a given Feynman diagram 
before the expansion in masses and true external momenta is made. The artificial external 
momenta have to be introduced in such a way that all spurious infrared divergences are re- 
moved, and the resulting Feynman integrals are calculable. Satisfying these two requirements 
is rather cumbersome in practical multiloop calculations. In addition, the condition that the 
IR divergences do not appear restricts considerably the power of the approach, since for com- 
plicated diagrams this requirement prevents one from reducing a given Feynman diagram to 
a simpler one. 

The latter problem was completely solved with elaborating a special technique of subtrac- 
tion of IR divergences — the .Reoperation [|]]. This method allows one to express (though 
in a rather involved way) the UV counterterm of every (h+l)-loop Feynman integral in 
terms of divergent and finite parts of some properly constructed h-loop massless propagators. 
Unfortunately, in practical applications, the use of the -Reoperation requires either many 
manipulations with individual diagrams or resolving a lot of non-trivial problem-dependent 
combinatorics (see, e.g. 0, §[). 

In our approach, the infrared rearrangement is performed by introducing an artificial 

mass rather than an artificial external momentum. A single mass parameter is added to 

2 In any meaningful renormalization prescription, counterterms are polynomial in external momenta, but 
not necessarily in masses. 
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each denominator of a propagator in each Feynman diagram. Consequently, no spurious IR 
divergences can appear. Next, an expansion in all the particle masses (except, of course, 
the auxiliary one) and external momenta is performed. The integrals one is left with have 
relatively simple form: They are completely massive tadpoles, i.e. Feynman integrals without 
external momenta and with only a single mass inserted in all the propagators. As a result, 
the problem of evaluating h-loop UV counterterms eventually reduces to a computation of 
divergent parts of h-loop completely massive tadpoles. 

At two loops, simple formulae for such Feynman integrals have been known since long 
ago [0]. However, no explicit formulae for three-loop massive tadpoles have been published 
so far. The available recursion algorithms M, [| based on the integration by parts method 



10 , |TTf are quite involved. 

The basic idea of our algorithm is to determine the pole part of a massive tadpole by 
expanding a properly chosen two-loop sub-integral with respect to its large external momen- 
tum being a loop momentum in the initial three-loop integral. Eventually, we have been 
able to construct relatively simple explicit formulae for all the necessary three-loop scalar 
integrals. 

The algorithm described in the present paper was used at the two loop fl2| and three 
loop |13| levels for calculating QCD anomalous dimensions of effective operators mediating 
B — > A s 7 decay. 

In principle, our method is applicable at the four-loop level, too. In this case, the problem 
eventually amounts to expanding a three-loop massive sub-integral of the propagator type 
with respect to its large external momentum. The algorithm for calculation of such three- 
loop integrals has been known since long ago (for a review see |14[]) and its computer algebra 
implementation has been recently achieved fl"5|| . 

Very recently, an alternative algorithm was developed by van Ritbergen, Vermaseren and 
Larin, and applied for evaluating four-loop contributions to the QCD (3 function and quark 
mass anomalous dimension [|l6j . Their approach amounts to using an identical to ours version 
of the IR rearrangement which reduces the calculation of UV renormalization constants to 
calculation of massive tadpoles. The difference appears at the stage of tadpole evaluation. 
The authors of [|16|] have succeeded in creating "special routines (...) to efficiently evaluate 4- 
loop massive bubble integrals up to pole parts in e and correspondingly of the 3-loop massive 
bubbles to finite parts." Eventually, all the diagrams have been reduced to two master ones. 
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Our paper is organized as follows: In the next section, we give general arguments which 
justify the use of an artificial mass parameter as an infrared regulator in all the propaga- 
tors, including propagators of massless gauge bosons. This is allowed so long as we are 
interested only in the UV-divergent parts of regularized Green's functions (with all UV 
su&divergences being pre-subtracted) . In section 3, we describe our algorithm for evaluating 
scalar integrals up to three loops. In section 4, we present some more details concerning 
calculation of nontrivial three-loop integrals. Section 5 contains two examples of relations 
between renormalization constants and (3 functions or anomalous dimensions up to three 
loops. Appendix A is devoted to reduction of tensor integrals to scalar ones. Appendix 
B summarizes expressions for "trivial" integrals, i.e. the ones which reduce to products of 
lower-loop integrals. Appendix C describes the expansion of one-loop self-energy integrals 
at large external momentum, which constitutes an essential element in calculating nontrivial 
three-loop integrals. Finally, appendix D contains a useful relation between tensor and scalar 
one-loop integrals in different numbers of dimensions. 

2. Decomposition of propagators 

The starting point of our procedure is a certain exact decomposition of propagators. For 
a scalar propagator belonging to a given Feynman diagram, it has the following form: 

1 _ 1 M 2 -p 2 -2qp-m 2 1 

(q + p) 2 — M 2 ~ q 2 -m 2 + q 2 - m 2 (q + p) 2 - M 2 ' ^ 

Here, p is a linear combination of external momenta in the considered diagram, q stands 
for a linear combination of loop momenta, and M denotes the mass of the particle. The 
artificial mass parameter m is introduced to regularize spurious infrared divergences. It is 
the same in all the propagators and all the diagrams. 

The contribution of the considered propagator to the overall degree of divergence of a 
diagram is Au = —2. The decomposition has been performed in such a way that the first 
simple term in the r.h.s. of eqn. ([]]) gives Au = —2, while the second, more complicated 
term gives Aui = —3. Moreover, the very last term in eqn. (0) has the same form as the 
original propagator. Thus, we can decompose it in an identical way. Doing so several times, 
we decompose the original propagator into a sum of terms with very simple denominators 
(depending only on loop momenta and the mass parameter m), and a more complicated 
term whose contribution to the overall degree of divergence is arbitrarily low negative. For 
instance, after three steps of decomposition, the exact expression for the original propagator 
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reads 



M 2 — p 2 — 2qp (M 2 - p 2 - 2qp) 2 



(q + p) 2 — M 2 q 2 — m 2 (q 2 — m 2 ) 2 (q 2 — m 2 ) 3 

m 2 m 4 — 2m 2 (M 2 — p 2 — 2qp) 

(q 2 — m 2 ) 2 (q 2 — m 2 ) 3 

(M 2 - p 2 - 2qp - m 2 f 



(2) 



(q 2 -m 2 f[(q + p) 2 - M 2 \ 
Here, the last term gives Au = —5 contribution to the overall degree of divergence of a 
diagram. 

In the following, we shall assume that the theory we consider is given by an (effective) 
lagrangian which does not contain non-negligible operators of arbitrarily high dimension, i.e. 
we assume that dimensionality of our operators is bounded from above. In such a case, any 
particular Green's function has a certain maximal degree of divergence. Consequently, we 
can always perform so many steps in the propagator decomposition, that the overall degree 
of divergence of any diagram in this Green's function would become negative if any of its 
propagators was replaced by the last term in the decomposition. We are then allowed to 
drop the last term in each propagator decomposition. It does not affect the UV-divergent 
part of the Green's function (after subtraction of sub divergences). 

It is important to note that each term in the propagator decomposition satisfies the 
criteria a full propagator should satisfy in the proof of Weinberg's theorem |17| |. This allows 



to apply degree-of-divergence arguments for diagrams where propagators are replaced by 
particular terms in their decomposition. 

A further simplification can be achieved by noticing that terms containing m 2 in the nu- 
merators (like the second line in the r.h.s. of eqn. (Q)) contribute only to such UV-divergent 
terms which are proportional to m 2 . These terms are local after subtraction of subdiver- 
gences. They must precisely cancel similar terms originating from integrals with no m 2 in 
propagator numerators. No dependence on m 2 can remain after performing the whole calcu- 
lation, because the propagators are decomposed exactly, i.e. they are actually independent of 
m 2 . This observation allows to avoid calculating integrals with m 2 in propagator numerators. 
Instead of calculating them, one can just replace them by local counterterms proportional to 
m 2 which cancel the corresponding (sub) divergences in integrals with no m 2 in propagator 
numerators. In effect, the practical calculation is made only with propagators replaced by 
such terms in the decomposition which contain no m 2 in the numerators (like the first line 
in the r.h.s. of eqn. (§)). Nevertheless, the final results for the divergent parts of Green's 



functions are precisely the same as if the full propagators were used (after subtraction of 
sub divergences) . 



In our earlier paper fL2| , we gave somewhat different arguments for using a single mass 
parameter as an infrared regulator in calculating (3 functions and anomalous dimensions. The 
present considerations might be more convincing, because the propagator decomposition we 
discuss here is exact. Thus, m? can be kept arbitrary all the time. One does not need to 
consider the m 2 — > limit and worry about its commutativity with Feynman integration. 

For particles with spin other than zero, the decomposition is applied only to denominators 
of their propagators, provided they are the same as in the scalar propagator. Our algorithm 
is not applicable in theories where kinetic terms differ from the canonical ones, as e.g. in the 
Heavy Quark Effective Theory ||18|| . 

As we have explained, one does not need to calculate Feynman integrals containing m 2 
in propagator numerators, so long as extra counterterms proportional to m 2 are introduced. 
Such counterterms may not preserve symmetries of the theory. Fortunately, the number of 
these counterterms is usually rather small, because their dimension must be at least twice 
smaller than the maximal dimension of operators in the considered (effective) lagrangian. 
For instance, the QCD lagrangian is built out of operators of dimension less or equal 4. 
There is only a single possible gauge- noninvariant counterterm of dimension 2. It reads 

l -Z x m 2 GlG a ^ (3) 

i.e. it looks like a "gluon mass" counterterm.^] At one loop, we find (using the Feynman- 
't Hooft gauge and the MS scheme in D = 4 — 2e dimensions) 

z * = -ils (iV + 2/ >- (4 > 

where N is the number of colors and / is the number of active flavors. This counterterm 
cancels gauge-noninvariant pieces of integrals with no m 2 in propagator numerators. 



After dropping the last term in the propagator decomposition, the Feynman integrands 
one is left with depend only polynomially on particle masses and external momenta. These 

quantities can be factorized out. It remains to calculate integrals depending only on loop 

3 The "ghost mass" counterterm does not arise in the Feynman-'t Hooft gauge, due to the structure of 
the ghost-gluon vertex. 
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momenta and the artificial mass parameter m 2 . At one loop, the generic integral reads 

f dD( l (r\ 

J (2n) D (q 2 - m 2 ) n ' 1 ' 

Integrals arising at more loops are slightly more complicated, because they involve several 
loop momenta. Nevertheless, reducing any such integral to scalar integrals can be easily 
performed by contracting it with various products of metric tensors and solving the resulting 



system of linear equations. We have written a Mathematica [jl|J code which performs such a 
reduction up to three loops, for an arbitrary number of free Lorentz indices. Some elements 
of this procedure are outlined in appendix A. 

After the reduction of tensor integrals is performed, one is left with relatively small number 
of scalar integrals to calculate. It is convenient to use the euclidean metric in discussing their 
evaluation. The euclidean integrals arising at one, two and three loops are respectively as 
follows 

1 



r(2) 

- I nin 2 n3 

r(3) 

nin2n3n4ri5n6 



m 



-D+2n -S- 



7T 2 



d D q- 



1 



m -3D+2Sn 1 ^-^ 



[q 2 + m 2 ] n ' 

J [q\ + m 2 ] n i[g| + m 2 ] n 2[(g 1 - q 2 ) 2 + m 2 ] n * 

J d D gi d D q 2 d D q 3 x 



(6) 
(7) 



x 



[q\ + m 2 ] ni \q\ + m 2 } n2 [q 2 + m 2 } n3 [(q 2 - q 3 ) 2 + m 2 ] n4 [(g 3 - gi) 2 + m 2 ] n s[(gi - g 2 ) 2 + m 2 ] n& ' 



The chosen normalization makes them dimensionless. The integrals can be represented by 
scalar vacuum diagrams displayed in fig. p] with propagators raised to arbitrary integer powers 
Hi- The algorithm for their evaluation is described in the next section. 






Figure 1: Graphical representation of the integrals given in eqns. 



3. The algorithm for evaluation of scalar integrals. 

In this section, we assume we are interested in evaluating three-loop /3-functions or anoma- 
lous dimensions. We use the MS scheme with D = 4 — 2e dimensions. We need to be able to 
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evaluate I$ n2mnm&n6 up to 0{\), up to 0(1) and /W up to 0(e). The latter integral 

is known exactly from textbooks 

I« = r < w -g) (9 ) 

r(n) ' 1 j 

The two-loop integral I^{i 27l3 is totally symmetric under permutations of its indices. It 
reduces to a product of one-loop integrals when at least one of the indices is nonpositive (see 
appendix B for explicit formulae). On the other hand, when all the indices are positive, it 



can be found from the following relations [^T 



I (2) - — l(3ni -DU (2) + nM m - 1 {2) 1 

1 (n 1 +l)n 2 n 3 Sn ' nin 2 n 3 ^ lb2 V 1 (ni-l)(n a +l)n 3 ni(n 2 +l)(n 3 -l)J 



+ n 3 [l[ n l_ 1)n2 ( n3+1) 4 1 ) ( n2 _ 1)(n3+1 )]} (10) 



and 



where 



/ni - (l- e )(l-2e) h 2^ + Y S ' ] + ° (e) ' 

4 /"T/3 



5 2 = = r rfxln(2sin-) ~ 0.2604341. (12) 

9V3 Jo 2 

The recursion relation (|T0|) holds for rij > 1. It can be derived with use of integration by 
parts. The sum of indices in the integral on its l.h.s. is bigger than the sum of indices in 
each of the integrals on its r.h.s.. Thus, the recursion can be programmed into a computer 
algebra code just as it stands. Two-loop integrals one usually encounters in practice are then 
found within a fraction of a second. 



Figure 2: Graphical representation of the integral 7 11U11 . It is equivalent to the last diagram 
in fig. H] with Hi = ... = n 6 = 1. 

Let us now turn to the three-loop integrals I^ n2ri3nin , 6n6 - Here, we are interested in 
calculating only UV-divergent parts of them. It is instructive to subsequently consider three 
cases: 

• A: All the indices ni, ...,hq are positive. 
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• B: At least one of the indices is equal to zero. 

• C: None of the indices vanishes, but some of them are negative. 

(3) 

Case A: From among all the three-loop integrals with six positive indices, only /xiini * s 
UV- divergent. Integrals with larger positive indices have negative degree of divergence and 
no subdivergences. 

In order to calculate the UV-divergent part of j{j} m we write this integral as follows (see 
fig.0): 

4iim = ™ D+2 ^ J d D q * jV^m 2 ), (13) 

where 

Jd q±d q 2 ^ + m 2 ][ql + m 2 ][(gi - q 2 ) 2 + m 2 ][(q - g x ) 2 + m 2 ][(q - q 2 ) 2 + m 2 ] ^ ^ 
The latter integral is just the usual two- loop contribution to the wave function renormaliza- 
tion in the " A0 3 " theory. We show it in fig. [3]. It is a finite diagram, because it has negative 
degree of divergence and no subdivergences. Consequently, a finite-volume integration in 
eqn. (|H|) cannot give a 1/e pole. Such a pole can only arise from integration over large q 2 
in eqn. ([131) . Therefore, knowing the behavior of Jinn(? 2 , m 2 ) at large q 2 is enough to find 



(3) 

the UV-divergent part of Iiiiin- 



Figure 3: Graphical representation of the integral Jnm(q 2 , m 2 ). 



The two-loop integral Jnm{q 2 , m 2 ) has the following expansion at large q 2 



AlUv 2 ,™ 2 ) 



m 

~ 2 



2\ 5-D 



q»m \ q 



m 



6C(3) + Ol — 1 +0{t 



Inserting this result into eqn. ( [13|) and introducing an infrared cutoff A one finds 



r(3) 



1 



'linn 



oo In 
d(q 2 " 1 q 



[q 2 + m? 



Jiiin(? 2 , m 2 ) + (finite terms) 



m\ & 2((3) 



A 
2C(3) 



+ (finite terms) 



(finite terms). 



(15) 



(16) 



The way we have found the UV-divergent part of /}iim shows the basic idea for calcu- 
lating all the nontrivial three-loop integrals in the cases B and C. The UV-divergent parts 
of these integrals can be found by choosing some two-loop subdiagrams of the last graph in 
fig. [1] and considering their behavior at large external momenta. If the considered two-loop 

(3) 

subdiagram is finite, the calculation proceeds analogously to the case of / 111U1 . If it is diver- 
gent, a subtraction of the UV divergence needs to be performed. We describe this in more 
detail below. 

Case B: Now, we consider the case when at least one of the indices of I^)„ n ^ Anc „„ is equal 
to zero. Without loss of generality, we can assume that the vanishing index is n 6 . This 
is because of the tetrahedron symmetry: The last diagram in fig. |1] has the topology of a 
tetrahedron. Symmetries of a tetrahedron can be described as certain permutations of its 
edges. Such permutations of the indices (n 1; ...,n 6 ) leave our integral invariant. 
When 77-6 = 0, the integral can be written as 

4 3 i 2 n 3 n4n 6 o = m~ D+2n3 tC% J d D qj-^—^G 2 (n u n 5 , ^LJ G 2 (n 2 , n 4 , ^ , (17) 

where the one-loop integral G 2 is given by 

2 ( 2 ' q 2 ) I ^[p 2 + m 2 ] kl [(p — q) 2 + m?\ k2 ^ ^ 

The diagram corresponding to eqn. (|T7|) is shown in fig. |J. 

ni n 2 



Figure 4: Graphical representation of the three-loop integral with one nonpositive index 

Calculation of the integral (|17D depends on what values are taken by indices of the one- 
loop sub-integrals denoted by G 2 .f\ One can distinguish three situations: 

• B.l: When at least one of these indices [n\ 1 n 2 ^nA t or n§) is nonpositive, then the 

three-loop integral reduces to a product of one- and two-loop tensor integrals. The 

4 However, it does not make any qualitative difference whether the index 77.3 in the final integration is 
positive or not. 
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latter integrals can be easily reduced to scalar integrals which we are already able to 
calculate. Final formulae for such three-loop integrals are given in appendix B. 

B.2: When nx,ri2, 714, > and both integrals G2 are convergent (i.e. rii+n^ > 2 and 
ri2 + n,4 > 2), the calculation proceeds analogously to the case of Inim- We expand 
the convergent integrals G2 at large q 2 



1 s 

m 



G2 ^ = E 



TYl 2 \ T ( Til 2 \ r 6 

— j a(h, k 2 ,r) + I —J b(k 1 ,k 2 ,r) 



(19) 



« ^ r=0 

Only a few lowest terms in this expansion affect the pole part of the considered 
three-loop integral. Explicit expressions for a(ki,k 2 ,r) and b(ki,k 2 ,r) are given in 
appendix C. 

B.3: When rii, 112, n.4, > but one or both integrals G2 are divergent (i.e. rii = = 
1 and/or n 2 = = 1), we need to split the integral G2 into its pole and convergent 
parts 

G 2 (ki, fea, "^pj = -JikMz + G r 2 en (h, h, pj ■ (20) 

Inside the three-loop integral, the pole part of G2 is multiplied by a two-loop integral. 
Thus, we already know how to calculate its contribution to the UV-divergence of the 
three-loop integral. On the other hand, the " renormalized" part of G2 can be treated 
analogously to the case B.2, i.e. in the same way the whole G2 was treated when it 
was convergent. Expansion of the "renormalized" integral G r 2 en at large q 2 is identical 



as in eqn. (jOj), except for that a(ki, k 2 ,r) is replaced by 



a ren (/ci, k 2 , r) = a(k x , k 2 , r) - -5 lkl 5 lk2 5 0r . (21) 

One should not naively expect that a ren (/c 1 , k 2 ,r) contain no poles in e. Actually, they 
do contain simple poles which cancel with the poles of b(k\, /c 2 , r) in the expression for 

Qren 



Case C: Now, we consider integrals with some negative indices but with none of them 
equal to zero. Using the tetrahedron symmetry, we can assume without loss of generality 
that one of the negative indices is n^. Then, we consider two distinct situations: 

• C.l: When any of the indices 711,712,^4 or n?, is negative, the three-loop integral 
reduces to products of one- and two-loop integrals, similarly to the case when n§ = 0. 
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The explicit formulae given in appendix B apply both when n§ vanishes and when it 
is negative. 

• C.2: When all the remaining indices are positive or the only other negative index 
is n 3 , we can still represent the considered three- loop integral by the diagram shown 
in fig. £|. However, instead of the scalar one- loop integrals G 2 , we encounter tensor 
one-loop integrals. This does not lead to any real difficulty, because we are able to 
reduce tensor integrals to scalar ones. Nevertheless, the amount of necessary algebra 
can be drastically reduced when one makes use of certain tensor identities discussed in 
appendix D. 

In the above considerations, we have described a complete algorithm for calculating pole 
parts of the integrals defined in eqns. However, obtaining final formulae for nontrivial 

three-loop integrals in the cases B.2, B.3 and C.2 requires discussing a few more subtle points 
in their evaluation. This is what the next section is devoted to. 



4. More on nontrivial three-loop integrals. 

Let us first derive our final expression for the three- loop integrals in the cases B.2 and 
B.3. In both these cases, the considered three- loop integral can be written as 

Cn 5 o = rn-^ tt-S / d D q [q2 _^ m2]n Gr (rn, n 5 , ^ G r 2 en (n 2 , n 4 , ^ 



1 (2) 1 (2) 1 



(1) 
3 



2\ / ™2\- n 3 00 / m 2\ »3+n+r 2 -3 



(fimte terms) + f ^y o + ^ 



x 



(TV? \ 6 / Ul^ \ ^ 6 

— J + b(n 1 ,n 5 ,r 1 )b(n 2 ,n i ,r 2 ) i—j 



+ [a ren (ni, n 5 , n)b(n 2 , n 4 , r 2 ) + 6(m, n 5 , r 1 )a ren (n 2 , n 4 , r 2 )\ 



( 



n\ 2e' 

m \ 



\Q 2 J 

Similarly to eqn. fllED , an arbitrary infrared cutoff A has been introduced here. Assuming 
that A 2 > m 2 , we can expand in eqn. (p2]) 

L , my" 3 g (-n3)(-n3-l)...(-n 3 -fc + l) g ^ ^ ^m»\ fc 

(23) 
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After performing trivial integrations, we arrive at the following result: 
4 3 i 2 n 3 n4n 5 o = ( fini te terms) + — ^- E 



ri-i 

V 2 / r 1 ,r 2 ,fc=0 





/ 9 


T ) ( 




) /m 2N 






)' + 


U 2 , 



X 



2\ 3^ 



&(wi,n 5 ,ri)6(w2,W4,r 2 ) / m_ 
n 3 - 2 + ri + r 2 + k + e \ A 2 J ' n 3 - 2 + r x + r 2 + k + 3e \ A 2 

^a ren (ni,n 5 ,r 1 )b(n 2} n4,r 2 ) + b(n 1} n 5} ri)a ren (n 2 ,n 4} r 2 ) ( m 2 \ 2e 
n 3 - 2 + n + r 2 + fc + 2e \ A 2 / 

+ -^lni^ln 5 4,2 ) n 4 n 3 + -^lna^^ingna _ ^ ^lrn ^ln 2 ^ln 4 ^ln 5 ^ ■ ( 24 ) 

The curly bracket in the above equation contains no 1/e poles unless n 3 — 2 +n +r 2 + = 0. 
Verifying this requires a short calculation, because a ren (ki, k 2 , r) and b(ki,k 2 ,r) do contain 
simple poles in e. Thus, for n 3 — 2 +r\ + r 2 + 7^ 0, one needs to expand the denominators to 
(9(e) and check that the potential 1/e contributions to J*- 3 -* "miraculously" sum up to zero, 
due to 

a Ten (k 1 ,k 2 ,r) + b(k 1 ,k 2 ,r) = 0(l). (25) 

(3) 

Our final expression for /„ in2n3 „ 4n5 o m the cases B.2 and B.3 is thus given by a finite sum 
(from now on we set A = m for simplicity^ ): 



J ?~o = (finite terms) + £ [ 



2—U3 2— ri3— ri 



X 



x |a ren (ni, n 5 , ri)a ren (n 2 , n 4 , r 2 ) + n 5 , ri)6(n 2 , n 4 , r 2 ) 

+ ^ [a ren (ni, 715, ri)6(n a , n 4 , r 2 ) + 6(n x , n 5 , ri)a ren (n 2 , n 4 , r 2 )] 

+ -5lni^ln5^n2 ) n 4 n3 + ~^n 2 ^lnJn^n s n 3 ~ ^In^lni^ln^lnJn} • (26) 

Let us now turn to the most complicated case C.2. In this case, the indices n%, n 2 , n 4 and 
n 5 are positive, while the index n 6 is negative. Using the tensor identities given in the end 
of appendix D, we express [(qi — q 2 ) 2 + m 2 ] - ™ 6 in terms of symmetric and traceless tensors 



-716 k [t/2] p fc-i+p / \ / \ / \ / . 

[( ?1 - ? 2) 2 +m 2 ]- = EEE E E T M 7 P 



X 



' ' ' ' iql + m 2 )- n ^- h {ql + m 2 ) k -^- l \ qi -q 2 f-^\ (27) 



4^p!(i-2p)!(i + 2-2p-e) /t 



5 One could keep A arbitrary and verify that the pole part of 1^ is independent of this parameter. This 
can serve as a useful cross-check against misprints in the explicit expressions for a(ki, k2,r) and b(ki, k2,r) 
in appendix C. 
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where (x) n denotes the Pochhammer symbol 

(x) n = x(x + l)(x + 2)...(x + n — 1) 
Consequently, we can write 



T(x + n) 
fix) 



(28) 



-n 6 k [i/2] p k-i+p / \ / \ / 

J(3 ) = yyyy y ~n 6 M P 

" 6<0 fc=oi=o P =oh=o h=o V fc / V 1 ) \ H 



k—i+p 
h 



p!(z-2p)!(z + 2-2p-e), 



(29) 



where 



r(3)(n) 

'nin2n3?i4n50 



-3D-2n+2En, 



3£j r 

7r _i 2- / cPgx cPg 2 ^ D <?3 x 

(?1 • g 2 ) (n) 



[g 2 + m 2 ] n i [g 2 + m 2 ] n 2 [g| + m 2 ] n 3 [(g 2 - gs) 2 + m 2 } n * [(g 3 - gi) 2 + m : 



2W 



(30) 



The above integral is a generalization of /^j^^go = Inin 2 n 3 n A n 5 o considered in the 
case B. When ni,n 2 ,n 4 or n 5 is nonpositive, I^n"n 3n4ns o * s equal to a linear combination 
of reducible integrals considered in the cases B.l and C.l. The explicit form of this linear 
combination is given in the end of appendix B. On the other hand, when rii, n 2 , 77,4 and 
are all positive, the calculation of InimnsmnsO proceeds analogously to the cases B.2 and 
B.3. However, instead of the scalar integrals G 2 , we encounter tensor one-loop integrals with 
totally symmetric and traceless tensors in their numerators. Such one-loop integrals are in 
one-to-one correspondence with scalar one-loop integrals in larger number of dimensions. 
The appropriate relation is given in appendix D. Using this relation, one finds the necessary 
generalization of eqn. ([17|) 



r (3)(n) 



n\ri2n 3 n 4 nh 



= (n 4 ) n (n 5 ) n m~ 



4+2e-2n+2n 3 „-2+e / j4-2e 



7T 



(n) 



[q 2 + m 2 ] ns 



x G 2 ni, n 5 + n. 



m 



G 2 n 2 ,n i + n. 



m 



D=4+2n-2e 



(31) 



D=4+2n-2e 



(32) 



where (g • q)^ can be expressed back in terms of g 2 

Similarly to eqn. (120|) , we split the higher-dimensional G 2 into its pole and convergent 
parts 



G 2 h, k 2 + n. 



m 



D=4+2n-2 e <n+l)\ \ q< 



D=4+2n-2e 

(33) 
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Next, we expand the convergent part at large q 2 , as in eqn. (|T^) 

^2 



/ 777 

G T 2 en lh,k 2 + n,— 



D=4+2n-2e 



E 



0^ — *OQ n 



^pj A ren (kx, k 2 + n,2 + n,r) + (^\ B(k 1 , k 2 + n,2 + n, r) 



(34) 



where 

A ren {k x , k 2 + n,2 + n,r) = A{k 1: k 2 + n,2 + n,r) - 1 ru 5 lkl 5 lk2 5 r (35) 
The coefficients A(k\, k 2 , u, r) and B{k\, k 2 , u, r) are given explicitly in appendix C. 



At this point, we are ready to write down the desired generalization of eqn. (|26|) 

4ii£ 3 n4n50 = (finite terms) + 



x 



nin2ri3?i4n50 

(n 4 ) w (n 5 ) n (2 - 2e) w 1 2 +" a+n-ns-n / 

2»(l-e) n eT(2-e) ^ ^ ^ 2+n-n 3 -n-r 2 

x {A ren (ni, n 5 + n, 2 + n, ri)A ren (n 2 , n 4 + n, 2 + n, r 2 ) 

+ \-B{ni,n b + n, 2 + n, r 1 )S(n 2 , n 4 + n, 2 + n, r 2 ) 
o 

+ [^en^ ni) ns + n? 2 + n, ri)-B(n 2 , n 4 + n, 2 + n, r 2 ) 
+S(n x , n 5 + n, 2 + n, r 1 )A ren (ra 2 , n 4 + n, 2 + n, r 2 )]} 



I [n/2] p n-2p j n-p-j / \ / n \ / . \ / 

E E E E E '" "~" / ' ' ^ ' "~ / '" / 



e 2«(n + 1) ^ ^ U t^o fo \ 1 J \ 3 J \ 1 



k 



X- 



p\{n-2p)\(n + l- p-e) p 



A A T {2) 
_i_x x r(2) 

-r«ln2 O ln 4 -'(n 1 -i-0(™5-i+0("-3- fc ) 



A A A A (2 - 2€) ra (2 - 6) ra r(n 3 -n-2 + e) 

4tel ^ Wlns ^fR ffa) • (36) 

The above equation is the main result of the present paper. It gives us pole parts of all 
the nontrivial scalar three-loop integrals I^ 3 \ i.e. those which do not reduce to products of 
lower-loop integrals. When n = 0, it reduces to eqn. 
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5. From renormalization constants to /5-functions and anomalous dimensions. 

In the preceding sections, we have described an algorithm for calculating pole parts of 
Feynman diagrams. Using our formulae, one can find all the MS-scheme renormalization 
constants in a given theory, up to three loops. In the present short section, we give two 
examples of relations between three-loop renormalization constants and beta functions or 
anomalous dimensions. 

Here, we depart from the MS scheme and assume that the renormalization constants 
(calculated in the framework of dimensional regularization) can contain arbitrary finite terms. 
However, we assume that these finite terms are renormalization-scale independent. 

For instance, let us consider renormalization of the gauge coupling g in some Yang-Mills 
theory 

g BARE = // Z 9 g, (37) 

where \i is the renormalization scale. The renormalization constant Z g has the following 
expansion in powers of the renormalized coupling g: 

( K U \ I K 12 K 22 \ ( K 13 K 23 K 33 \ 

Z9 = 1+92 ( kM + ^) + ^(^ + ^ + ^) +9 »( k » + ^ + ^ + ^ + ... (38) 

Some coefficients in this expansion are given in terms of the others, which follows from 
locality of UV-divergences 

^ = \{^f (39) 

K 23 = H^V 2 _ I k oi (k ii )2 
3 2 K ' 

From scale-independence of g BARE one can derive the following expression for the (3- 
function in terms of k % ^\ 



9 7 +- 



P( g ) =n^- = 2K U g 3 +UK 12 - 12« 01 « n ] g 5 +[6K 13 - 22k 1 V 2 - 22k 01 k 12 + 54k u (k 01 ) 2 

(40) 

As another example, let us discuss the anomalous dimension matrix of a set of (possibly 
dimensionful) couplings Cj which linearly mix under renormalization 

C»*** = Y,C i Z ij =(p T z).. (41) 
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Let us assume, that the renormalization constant matrix Z depends on a single gauge cou- 
pling g. Then it reads 

z = i + ! ? (a- + *pj +9 * (a« + v + ^ ) ( a ° 3 + T + ? + + • • (42) 

Some coefficients in the above expansion are given in terms of the others, which follows 
from locality of UV- divergences 

a 22 = ^(a 11 ) 2 + K n a n 

a 33 = I(a 11 ) 3 + K 11 (a 11 ) 2 + ^( K 11 )V 1 

a 23 = ia n a 12 + -a 12 a n - V(a n ) 2 - -a n a 01 a n - i/c^o" + -/^a 12 + f -k 12 - -k 01 k u 
3 3 6 V 7 3 3 3 V3 3 

Scale-independence of (J BARE implies that the renormalized couplings Cj satisfy the fol- 
lowing renormalization group equations 

= 7 T C, (44) 
where the anomalous dimension matrix 7 has the following expansion in powers of g 



7 = la g + g 4a — za a — za a — 4k a — 4k a +5 



-12 oaOI-11 oaIIaOI /i ,.11 -01 „ ,01 A ll 



,6 



6a 13 -4a 12 a 01 -2a 01 a 12 



-4a 02 a n - 2a n a 02 + 2a 01 a 11 a 01 + 2a n (a 01 ) 2 + 2(a 01 ) 2 a n + 4 K n (a 01 ) 2 - 8K n a 02 
+4/ t ° 1 a 11 a 01 + 4K 01 a 01 a n - 8K 01 a 12 - 8K 12 a 01 - 8K 02 a n + 24A 11 a 01 + 12( K 01 ) 2 a n ] + ... 

In the MS scheme, equations (|40|) and ( f£5|) become much simpler 

/%) = 2/t 11 ^ 3 + 4/t 1 V + 6K 1 V + ... (46) 
7 = 2a 1 V + 4a 1 V + 6a 1 V + ... (47) 

However, using the pure MS scheme may not be possible in some effective theories where 
so-called "evanescent operators" arise in dimensional regularization. This is why the more 
general relations (f4"0|) and (|4"5|) have been presented here. 
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6. Summary. 

We have described an algorithm for calculating UV-divergent parts of arbitrary Feynman 
diagrams. A common mass parameter has been used to perform the infrared rearrangement. 
Explicit formulae for all the necessary scalar integrals up to three loops have been given. 

The main idea in calculating nontrivial three-loop integrals was considering some of their 
two-loop subintegrals and expanding them at large external momenta. In the end, some 
details have been given on relations between UV- divergences and /3-functions or anomalous 
dimensions. 
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Appendix A 

This appendix is devoted to reduction of tensor integrals to scalar ones. We are interested 
in Feynman integrands depending only on loop momenta and the artificial mass parameter 
m 2 . The integrals arising at one, two and three loops have the following form (in the euclidean 
metric): 

1 " JW "^ J [p2 + m 2]n > ^ 

(T (2)ki \ = 2D -k-l+2j2ni _-£> f d D pd D q Pn 1 ...Pn k q vl ...q Vl 

(T (3)Mm \ _ ^3D-k-l-m+2j2n t v 

y- 1 n 1 n 2 n i n i n 5 n e ,) m...li k vi...v l p 1 ...p m ~ " b " A 

d D p d D q d D r p^.-.p^ q vl -q Vl r pi ...r Pm 

[p 2 + m 2 ]™ 1 [q 2 + m 2 } 712 [r 2 + m 2 } n * [(q - r) 2 + m 2 ]™ 4 [(r - p) 2 + m 2 ] n5 [{p - q) 2 + m 2 ] n e ' 

(50) 
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Such integrals are proportional to linear combinations of products of metric tensors. For 
instance, 

(-^7iiTi2"3 )p u P cr -^1 9pu9pa -^2 (fiVpfiW + 9^9u P )- (51) 

The tensors g^, p g U a and g^Qvp are multiplied by the same coefficient F 2 in the above equation, 
due to an obvious symmetry. In a computer algebra code, such a symmetry can be verified 



by checking that contractions of the l.h.s of eqn. (|5l|) with g w g va and g^g^p are identical. 

The equations for the coefficients F\ and F 2 are found by contracting the tensor integral 
with g^gpa and g w g va 

f D 2 F 1 + 2DF 2 = X, 

\ DF 1 + D(D + 1)F 2 = X 2 K ° > 



where 



- 2D -4+2j> -d f d D p d D q p 2 q 2 



J [p2 +m 2]ni[ g 2 + m 2]n 2 [(p_ g )2 +m 2jn 3 ' ^ ; 

-2d-4+2J> -d /" rf> d D q (p ■ qf 



' [p 2 + m 2 ] n i[g 2 + m 2 ] n2 [(p - g) 2 + m 2 ] n 3 v ' 



Consequently, 



F x \ ( D 2 2D \ 1 f X 1 



\f 2 j ~ \ d d{d + 1) ; ^ x 2 ; ' (55) 

The above matrix inversion is most easily done perturbatively in e, after substituting D = 4 — 
2e. This makes the computer program much faster, which is important for more complicated 
tensor integrals where larger matrices need to be inverted. 

The integrals like Xi and X 2 are easily reduced to the standard scalar integrals (0), with 
help of the identities 

p 2 = (p 2 + m 2 )—m 2 , (56) 
q 2 = ( q 2 + m 2 ) - m 2 , (57) 
p-q = -{(p 2 + m 2 ) + (q 2 + m 2 ) — [(p — q) 2 + m 2 ] — m 2 }. (58) 

Specific values for the indices rii can be substituted only after all these operations are per- 
formed. 



Appendix B 

In this appendix, we give explicit formulae for the trivial integrals, i.e. for the two-loop 
integrals which reduce to products of one-loop ones, and for the three-loop integrals which 
reduce to products of one- and two-loop ones. 
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When at least one of the indices of the two-loop integral I$ n2n3 is nonpositive, the integral 
reduces to a product of tensor one-loop integrals. However, a simple expression for such an 
integral can be also obtained from a general formula for a two-loop integral with one massless 
and two massive lines [JJ 

d D qi d D q 2 



r(2) = _ V ( ~™ 3 ) TO -2o+2(n 1+ n 2 +fc) -D f 

V» 2 n 3 jL[ k j m J 



[qf + m 2 ]^ [ql + m 2 \ n * \{q x - q 2 ) 2 } 



21 —k 



^ ( -n 3 \ r(f + fc)r(m - fc - f )i> 2 - k - f )r(m + n 2 - k - d) 



In the above equation, we have assumed that the nonpositive index is n%. This could have 
been done without loss of generality, because I$ n2ri3 is totally symmetric under permutations 
of its indices. Equation (|59"D implies that I^ n2n3 vanishes when more than one of its indices 
is nonpositive. 

Let us now turn to the three-loop integrals ^n 2n3n4 „ 5n6 considered in the cases B.l and 
C.I in the main text. They have nonpositive index uq and, in addition, there is one more 
nonpositive index among ni,n 2 , and 715. Symmetries of the diagram shown in fig. f| allow 
to assume without loss of generality that the other nonpositive index is 715. The remaining 
indices can be arbitrary integers. In such a case, the three loop integral is expressible in 
terms of tensor one- and two-loop integrals as follows: 

-m -ng 35 h —ns—js —ne—je / \ / \/.\/.\ 

r(3) = V V V V V V (— 1V' 5+J ' 6 2 i5+ie 

nin2n:jri4,n5n(j / , / j / ^ / j / / / , V V " 
J5=0 j 6 =0 i 5 =0 i 6 =0 fc 5 =0 fc 6 =0 



~n 5 \ I -716 \ I 35 36 
35 I \ 36 I \ i5 I \ 16 



X 



/ -ns-js \ / -ne-ie \ / r T (l)(i 5 +i 6 ) > \' il -' ii 5+<6 ( T (2)i 6 i 5 \ / fifj x 

I k 5 J \ fc 6 i V 1 n 1 -k 5 -k 6 J \ J (n 2 +n 6 +i6+&6)(n3+n5+is+fc 6 )n4; Ml .,^ is+i6 ^ UU J 

The tensor integrals appearing in the above equation have been defined in appendix A. 

In the end, we consider the integral IninXsmnsO defined in eqn. ( p0|) in the reducible case, 
i.e. when ni,n 2 , n.4 or n 5 is nonpositive. In such a case, we can calculate ^nansrunso ^ 
expressing it as a linear combination of the integrals considered in the previous paragraph 



[n/2] n-2p k p k-i+p / \ / \ / \ / 

j(3)(n) [ n-2p \ I k \ I P \ I k-i+p 



X 



k f\i/\h/\ h 



p=0 k=0 i=0 h=0 l 2 =0 

C_1 \p+n+i+l 1 +h r) \ 

x I L > a - j( 3 ) (Ql) 



In some of the integrals I® on the r.h.s. of the above equation, we may need to permute the 
first five indices using the tetrahedron symmetry. The fifth index must become nonpositive 
before eqn. (|60|) is applied. 



19 



Appendix C 

Here, we give explicit formulae for the coefficients a(ki,k2,r) and b(ki,k 2 ,r) in the ex- 
pansion (|T9| ) of G 2 at large q 2 . In section 4, we also need their generalizations to D = 2uo — 2e 
dimensional space, with arbitrary u. Thus, we write 



a(k l7 k 2 ,r) = A(k 1 , k 2 , 00 = 2, r) 
b(ki,k 2 ,r) = B(k 1 , h 2 , u) = 2, r). 



(62) 
(63) 



The above quantities are symmetric with respect to their first two arguments. More- 
over, we are only interested in positive k\ and k 2 in our application to three-loop integrals. 
Consequently, knowing A and B for 1 < k\ < k 2 is everything we need here. 

Using eqns. (18) and (A.l) of ref. ]22]]|] one finds (for 1 < k\ < k 2 ) 



A(k l ,k 2 ,u,r) 







(-l) r - fc i(r-l)!(fc 2 -r-l)!r(fci+fc 2 -r--4J+6) 
(fei-l)!(fc 2 -l)!(r-fci)!(fci+fc 2 -2r-l)! 

2(r-l)!(2r-fci-A; 2 )!r(fci+fc 2 -r-a;+e) 
(fci-l)!(A: 2 -l)!(r-fc 1 )!(r-A: 2 )! 



when r < k\ or ^ 2 < r < k 2 
when ki <r < kl + k2 
when r >k 2 



B(ki,k 2 ,u,r) 







(_!)■— fcl-fc2+"r(fci-r-e)r(fc 2 -r-e)r(r+e) 
(fci-l)!(fea-l)!(r-fei-Jfc 3 +w)!r(A!i+fea-2r-2e) 



(64) 
(65) 



when r < fci + — co> 
when r > ki + k 2 — u 

The coefficients A and S often contain simple poles in e. For convergent integrals, these 
poles are usually also present, but they cancel out in the final expression for G 2 . Thus, even 
when both G 2 in eqn. (|T7|) are convergent, one needs to carefully keep track of the 0(e) parts 
in A and B. 



Appendix D 

Here, we present a useful relation between tensor and scalar one-loop integrals in different 
numbers of dimensions. Let p( ai - a «) denote the only symmetric and tracelessj] rank n tensor 
which can be formed from a D- vector p (see eg. |23[|) 

[ n / 2 l I 1 \p I 

Here, S stands for the operator which symmetrizes with respect to all n indices and multiplies 
the result by l/n\. 

6 Elementary identities T(z)T(z + 1/2) = 2 1 ~ 2z r(2z) and T(z)T(l - z) = tt/ sin(Trz) have been used to 
simplify these expressions. 

7 Tracelessness of a rank n tensor means vanishing of any of its contractions with g^ v . 
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The tensors p( a i- a ») occur in a useful relation between one-loop tensor and scalar integrals 



in different numbers of dimensions |24 

m ~D+2k 1+ 2k 2 n -D/2 f d D p 



P 



(ai...a n ) 



[p 2 + m 2 ] kl [(p — q) 2 + m 2 ] k2 
{k 2 ) n q iai --- a - ) m- D - 2n+2kl+2k ^- D / 2 - n J d D+ 



V 



[p 2 + m 2 ] kl [(p — q) 2 + m 



21fc 2 +" 



(67) 



D=4+2n-2e 

The tensor one-loop integrals one finds in calculating scalar three-loop integrals are not 
given in terms of symmetric and traceless tensors. Instead, one encounters powers of scalar 
products of various D-vectors. However, these latter objects can be reversibly related to 
contractions of symmetric and traceless tensors. Let us define 



(p-9) ( "W ai --- a "W..««)- 



(68) 



In eqns. (A. 10) and (A. 15) of ref. ||11|| , the following relations have been given:^ 

(-l)o 



(p.g)W 

(p ■ <D n 



[n/2 

E 

p- 

[n/2] 



n\ 



E 



" (n+ 1 - p-e) p APp\(n-2p)\ 
1 n! 



{p 2 q 2 Y{p-q) n - 2p (69) 
{p 2 q 2 ) p {p-q) {n ~ 2p) (70) 



p= - (n + 2-2p-e) p 4v!(n-2p)! 

This appendix summarizes all the information on symmetric and traceless tensors necessary 
for deriving our final expressions (|29|) and (|36|) for the three-loop integrals in the most 
complex case C.2. 



8 There is a misprint in eqn. (A. 15) of ref. JliJ which we correct here. 



21 



References 

[1] J.C. Collins, Nucl. Phys. B92 (1975) 477. 

[2] A.A. Vladimirov, Teor. Mat. Fiz. 43 (1980) 210. 

[3] K.G. Chetyrkin, A.L. Kataev and F.V. Tkachov, Nucl. Phys. B174 (1980) 345. 

[4] K.G. Chetyrkin and V.A. Smirnov, Phys. Lett. B144 (1984) 419. 

[5] K.G. Chetyrkin, Phys. Lett. B390 (1997) 309. 

[6] K.G. Chetyrkin, Phys. Lett. B391 (1997) 402. 

[7] J. van der Bij and M. Veltman, Nucl. Phys. B231 (1984) 205. 

D.J. Broadhurst, Z. Phys. C54 (1992) 599. 
[9] L.V. Avdeev, Comput. Phys. Commun. 98 (1996) 15. 
[10] F.V. Tkachov, Phys. Lett. B100 (1981) 65. 

[11] K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. B192 (1981) 159. 

[12] M. Misiak and M. Munz, Phys. Lett. B344 (1995) 308. 

[13] K. Chetyrkin, M. Misiak and M. Munz, Phys. Lett. B400 (1997) 206. 

[14] V.A. Smirnov, Mod. Phys. Lett. A 10 (1995) 1485. 

[15] K.G. Chetyrkin, R. Harlander, J.H. Kuhn and M. Steinhauser, 
Nucl. Instrum. Meth. A 389 (1997) 352; 
preprint MPI/PhT/97-012, TTP97-11 (|hep-ph/ 9 7042221) . 



[16] T. van Ritbergen, J.A.M. Vermaseren and S.A. Larin, Phys. Lett. B400 (1997) 379; 
Phys. Lett. B405 (1997) 327. 

[17] S. Weinberg, Phys. Rev. 118 (1960) 239. 

[18] M. Neubert, Phys. Rept. 245 (1994) 259. 

[19] S. Wolfram, Mathematica, Addison- Wesley Publ. Comp., N.Y. 1988. 

[20] J. Collins, Renormalization, Cambridge Univ. Press, Cambridge, 1984. 

[21] A.I. Davydychev and J.B. Tausk, Nucl. Phys. B397 (1993) 123. 

[22] E.E. Boos and A.I. Davydychev, Theor. Math. Phys. 89 (1992) 1052. 

[23] A.V. Kotikov, Phys. Lett. B376 (1996) 240. 

[24] A.I. Davydychev, Phys. Lett. B263 (1991) 107; 
O.V. Tarasov, Phys. Rev. D54 (1996) 6479. 



22 



